Evolution of magnetism in Yb(Rhi_ x Co x ) 2 Si2 
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We present a study of the evolution of magnetism from the quantum critical system YbRh2Si2 
to the stable trivalent Yb system YbCo2Si2- Single crystals of Yb(Rhi_ I Co a; )2Si2 were grown 
for < x < 1 and studied by means of magnetic susceptibility, electrical resistivity, and specific 
heat measurements, as well as photoemission spectroscopy. The results evidence a complex mag- 
netic phase diagram, with a non-monotonic evolution of Tjv and two successive transitions for some 
compositions resulting in two tricritical points. The strong similarity with the phase diagram of 
YbRli2Si2 under pressure indicates that Co substitution basically corresponds to the application 
of positive chemical pressure. Analysis of the data proves a strong reduction of the Kondo tem- 
perature Tk with increasing Co content, Tk becoming smaller than TV for x ~ 0.5, implying a 
strong localization of the 4/ electrons. Furthermore, low-temperature susceptibility data confirm a 
competition between ferromagnetic and antiferromagnetic exchange. The series Yb(Rhi_ I Co I )2Si2 
provides an excellent experimental opportunity to gain a deeper understanding of the magnetism at 
the quantum critical point in the vicinity of YbRli2Si2 where the antiferromagnetic phase disappears 
(liv^O). 

PACS numbers: 71.10.Hf 71.27.+a 75.20.Hr 75.30.-m 



I. INTRODUCTION 

In recent years the heavy fermion (HF) metal 
YbRli2Si2 has been intensively investigated due to 
its proximity to an antiferromagnetic (AFM) field- 
induced quantum critical point (QCP) 1 . This proxim- 
ity leads to extraordinary non-Fermi-liquid (NFL) be- 
havior, such as the divergence of the electronic Sommer- 
feld coefficient^— the linear-in-T resistivity^ as well as 
diverging Griineisen ratios^ The reason of the NFL be- 
havior is assumed to result from quantum fluctuations 
with a strong ferromagnetic (FM) component^. More- 
over, a jump of the Fermi volume has been observed in 
Hall-effect^ii^ and thermopower— measurements, which 
could not be explained by the conventional 3D spin- 
density-wave (SDW) scenario^ 2 - The recent discoveries of 
an additional energy scale vanishing at the QCP, which 
does neither correspond to the Neel temperature nor to 
the upper boundary of the Fermi-liquid region 13 , and a 
large critical exponent a = 0.38 at the AFM phase transi- 
tion observed in low-temperature specific-heat measure- 
ments on a single crystal of superior quality^ have once 
again boosted the interest in YbRli2Si2- Theoretically, 
the unconventional quantum criticality in this material 
may be ascribed to a breakdown of the Kondo screening 
at the field induced AFM QCP^r-ii 

Besides strong experimental efforts, the determination 
of the magnetic structure using neutron scattering ex- 
periments was not yet successful due to the small size 
of the ordered moment (» 10~ 3 hb)^ and the very low 
Neel temperature Tjy = 72 mK. One way to overcome 
these difficulties would be the application of hydrostatic 
pressure, which results in a stabilization of the AFM or- 



dering and, therefore, in an increase of the ordered mo- 
ment and In addition, the appearance of a sec- 
ond phase transition, labeled XL, has been reported at 
p > 1 GPa.— An alternative to hydrostatic pressure is 
isoelectronic doping with cobalt which leads to chemi- 
cal pressure because of a smaller unit-cell volume, V. 
Hodges already showed earlier that YbCo2Si2 orders an- 
tiferromagnctically at TV. = 1.7 K with a sizable or- 
dered Yb moment of 1.4 /is.— Later on, susceptibility 
and inelastic neutron scattering experiments confirmed a 
stable trivalent Yb state with well defined crystal elec- 
tric field levelsjSi Several results on RCo2Si2 compounds 
proved Co to be non-magnetic, because of a strong Co- 
Si hybridization^ First experiments with small amounts 
of Co doping led already to very promising resultS) 23 i 24 
when compared to the equivalent hydrostatic pressure ex- 
periments on YbRh2Si2 J 19 ' 25 ' 26 We therefore studied the 
complete doping series Yb(Rhi_. E Co a ;)2Si2 in order to get 
a reliable picture of the evolution and the characteristics 
of the magnetic order. 

The paper is organized in the following way. In the 
first section the crystallographic analysis of all doped 
samples are given, followed by the presentation of the 
magnetic susceptibility results. The subsequent sections 
are dedicated to the electrical-resistivity data followed by 
the presentation of the specific-heat results. Finally, the 
experimental results of the photoemission spectroscopy 
are shown. From these data we derive the magnetic 
phase diagram of the complete series Yb(Rhi_ a; Co a ;)2Si2. 
Its close resemblance to the magnetic phase diagram of 
YbRli2Si2 under pressure confirms that the present alloy 
system provides an excellent access to a deeper under- 
standing of the unusual properties of YbRli2Si2. 
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FIG. 1: (Color online) Lattice parameters a and c of the 
Yb(Rhi_ :E Co a ;)2Si2 series. The slope of the data indicates 
positive chemical pressure with increasing x. Left inset: Pic- 
ture of typical single crystals for all Co concentrations. Right 
inset: c/a ratio versus unit cell volume. The lines represents 
linear fits to the experimental data. 

II. EXPERIMENTAL RESULTS 
A. Crystallographic parameters 

Single crystals of Yb(Rh 1 _ :c Co :; .)2Si2 were grown us- 
ing the indium flux-growth technique developed for 
YbRli2Si2 27 . Surprisingly, the growth parameters opti- 
mized for the pure Rh-based compound could be retained 
for the whole alloy series, enabling us to grow large single 
crystals up to pure YbCo2Si2 (see left inset, Fig. [1]). All 
samples were prepared in A^Oa-crucibles, which were 
placed in a second crucible of tantalum closed under ar- 
gon atmosphere using arc melting. The elements were 
heated up to T = 1500° C, followed by a slow cooling 
to T = 1000° C. After the growth, the single crystals 
were retrieved from the In-flux using hydrochloric acid for 
chemical etching. None of the crystals showed any indi- 
cation of surface destruction due to HC1 attack. Only for 
the Co concentration x = 0.27 the surface of the crystals 
appeared to be a little different, probably due to residue 
elements after the etching process (see left inset, Fig. [I}. 
However, after polishing the surface, different physical 
measurements led to the conclusion that the bulk prop- 
erties were not affected. 

To find out the crystal orientation, the single crys- 
tals were analyzed using Laue backscattering. The Laue- 
reflection pictures confirmed that the surface of the crys- 
tal is always perpendicular to the crystallographic c-axis. 
Therefore, the crystal growth in the a6-direction is fa- 
vored with respect to the c-direction. Subsequently, X- 
ray powder-diffraction patterns on crushed single crys- 
tals showed that the complete doping series is single 
phase with the expected ThCr2Si2-structure (space group 
I4/mmm), in agreement with the crystal structure of 
YbRli2Si2 and YbCo2Si2. 28 Since one aim of this work 



TABLE I: Lattice parameters for the complete series obtained 
from X-ray analysis and least square root fitting procedures. 



X 


a (A) 


c (A) 


c/a 


v (A 3 ) 


±0.01 


±0.001 


±0.002 


±0.001 


±0.1 


0.0 


4.011 


9.861 


2.458 


158.6 


0.03 


4.008 


9.855 


2.459 


158.3 


0.07 


4.003 


9.846 


2.460 


157.8 


0.12 


3.997 


9.836 


2.461 


157.1 


0.18 


3.988 


9.821 


2.463 


156.2 


0.195 


3.986 


9.820 


2.464 


156.0 


0.215 


3.984 


9.816 


2.464 


155.8 


0.27 


3.979 


9.804 


2.464 


155.2 


0.38 


3.961 


9.787 


2.471 


153.6 


0.58 


3.934 


9.764 


2.482 


151.1 


0.68 


3.911 


9.743 


2.491 


149.0 


0.78 


3.893 


9.729 


2.499 


147.5 


1.0 


3.852 


9.689 


2.515 


143.8 



is to compare the effects of hydrostatic and chemical 
pressure, the cell parameters a and c should both de- 
crease with increasing doping. Otherwise, the compari- 
son would be more difficult due to a non-isotropic change 
of the unit cell. Fortunately, for the complete series 
Yb(Rhi_ x Co 2; )2Si2 both lattice parameters a and c con- 
tinuously decrease (Fig. [1] and Tab. [I]). The ratio c/a 
shows only a weak linear increase with x with differ- 
ent slopes below and above x ps 0.38, the former being 
very small. Since hydrostatic pressure experiments evi- 
dence a constant c/a up to highest investigated pressure 
of p = 21 GPa,^ this slight change of c/ a has to be taken 
into account in the comparison of hydrostatic and chem- 
ical pressure results. 

The Co content was determined using energy disper- 
sive X-ray spectra performed on a scanning electron mi- 
croscope (Philips XL30) with a Si(Li)-X-ray detector. 
Three crystals of each batch with a polished surface were 
analyzed at several positions. The respective mean val- 
ues result in the given percentages with a relative error 
of less than 1 at.%. These real Co concentrations, which 
are used for x throughout this paper, are somewhat larger 
than the nominal concentrations due to different solubil- 
ities of the elements in the In-flux. 



B. Magnetic susceptibility 

In this section, we discuss the magnetic behavior of 
the Yb(Rhi_ a; Co a ;)2Si2 single crystals. All measurements 
were carried out using a Magnetic Property Measurement 
System (MPMS) of Quantum Design. We also present 
results of the dc-susceptibility down to T — 0.5 K, which 
where realized with a 3 He-option for the MPMS designed 
by the iQuantum Corporation. 

Dc-susceptibility (%) measurements on YbRli2Si2 re- 
vealed Curie- Weiss-law behavior between 200 < T < 
300 K with a large magnetic anisotropyi 2 ^ At high tem- 
peratures, the effective moments along a and c, /ig ff and 
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FIG. 2: (Color online) Temperature dependence of the mag- 
netic susceptibility for x = 0.18 and 0.27. The AFM tran- 
sition temperatures Tl and Tn are marked by arrows. The 
inverse susceptibility is exemplarily shown in the inset for 
x = 0.215 and 0.68 presenting paramagnetic Curie- Weiss be- 
havior above 150 K with /j, a g = 4.6/^b- 



/Ltg ff , equal the value of 4.54 /is corresponding to the J = 
7/2 multiplet of the free Yb 3+ ion. Because the strong 
easy-plane anisotropy is retained throughout the whole 
alloy, we focus here on the in-plane susceptibility results. 
In the temperature range 200 < T < 300 K, x(T) shows 
Curie-Weiss (CW)-like behavior with an effective mo- 
ment Hcf? — 4-6 ± O.I/Ub for all x (see Tab. Inland inset of 
Fig.EJ, as expected for a trivalent Yb state, and thus sup- 
porting the absence of a Co moment. The difference of 



the Weiss temperatures, 6^' 200 """ K 



©W ( see 



Tab. Ill]) , changes only slightly throughout the complete 
series, indicating a rather smooth variation of the leading 
crystal electric field (CEF)-parameter i? 2 (Stevens co- 
efficient) when going from YbRh 2 Si 2 to YbCo 2 Si 2 <2ML 
This is supported by our analysis of the CEF-scheme 
for x = l, 32 resulting in a T7 ground state similar to 
YbRb^Sia^ 

The first low-temperature ac-susceptibility measure- 
ments for x < 0.12 were performed and reported by 
Westerkamp et al<22 They already confirm the stabiliza- 
tion of the AFM ordering with increasing x. In addition 
to the increasing Tn, a second magnetic feature appears 
in the magnetically ordered region for x > 0.07. Vari- 
ous measurements on YbRh 2 Si 2 under hydrostatic pres- 
sure showed the appearance of a second anomaly, Tl, at 
p > 1 GPa . 19 i 26 Focusing now on the next higher doping 
x = 0.18, the two mentioned anomalies are clearly visible 
in Fig. [2] for a magnetic field B = 50 mT, applied per- 
pendicular to the hard magnetic c-axis. At Tn, a small 
change of slope in the temperature dependence marks the 
onset of AFM order; whereas, below Tl, the typical de- 
crease of x(T) is observed. The fact that the ordering is 
of AFM type can be deduced from the temperature and 
field dependence of the magnetic susceptibility, as well as 
from the shifting of Tn and Tl towards lower tempera- 



FIG. 3: (Color online) Low-temperature Weiss temperature, 
W ' ~ , for the complete series showing the increasing FM 
component of the exchange interactions with increasing x (left 
axis, circles). The effective moments in the same temperature 
region, pi^' 2_4K , indicate a decreasing Kondo screening for 
higher cobalt contents (right axis, diamonds). Dashed lines 
are guide to the eyes. The inset shows exemplarily the Curie- 
Weiss behavior together with the linear fit for two different 
concentrations. 



tures with increasing external magnetic fields (cf. Fig. [2]). 
For x — 0.18, an in-plane field of B = 0.3 T is sufficient 
to suppress the phase transition at Tl, while the mag- 
netic ordering has completely disappeared at B — 0.6 T. 
Increasing the doping to x — 0.27, Tl is not observable 
anymore for T > 0.5 K, whereas Tn is shifted towards 
higher temperatures. The T-dependence of x is still of 
AFM type but changes profoundly compared to the one 
observed for x = 0.18. Recent resistivity measurements 
on the x = 0.27 sample showed that no second phase 
transition appears down to T = 20 mK.— 

The important role of the FM fluctuations in YbRh 2 Si 2 
has been already mentioned. As pointed out by Ishida 
et alZ for YbRh 2 Si 2 and by Gegenwart et al£- for 
YbRh 2 (Sio.95Geo.o5), the quantum critical fluctuations 
have a very strong FM component. Among all other 
known quantum critical HF compounds it is also this fea- 
ture which makes the HF system YbRh 2 Si 2 unique. In 
the following, we discuss the evolution of the FM fluc- 
tuations above Tn upon increasing the chemical pres- 
sure. In systems with a stable Yb 3+ state, because 
of the small de-Gennes factor, magnetic inter-site ex- 
changes are much weaker (< 10 K) than crystal field ef- 
fects (< 100 K). Then, the Weiss temperature obtained 
from a fit to the low-temperature susceptibility data re- 
flects the sum of all exchange interactions (while the fit to 
the high-temperature data reflects the CEF effect) and is 
a good indicator for the dominant exchange interactions. 
We therefore fitted the measured inverse susceptibility 
at B = 0.5 T with a CW law between 2 < T < 4K, 
just above Tn (see inset of Fig. [3]) for the complete se- 
ries. The values of 6^' 2 4K and /i^' 2 ~ 4K deduced from 
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TABLE II: Effective moment, jU e ff, Weiss temperature, 0w, from the high-temperature (200 < T < 300 K) and the low- 
temperature (2 < T < 4K) CW fits for both the magnetic easy plane (ab) and the hard axis (c). The magnetic phase transition 
temperatures, TL and Tpj, were obtained from heat-capacity measurements or from the literature; * denote first-order phase 
transitions. 



X 


ab:200 — 300K / \ 
Meff (PS) 


,~.ab;200 — 300K i-tr\ 


e ^200-300K (K) 


ab;2 — 4K / \ 
Mcff (MS) 


e^ ;2 - 4K (K) 


Tl(K) 


T N (K) 




±0.1 


±10 


±10 


±0.2 


±0.2 


±0.03 


±0.03 


0.0 


4.5 


-9 


-215 


2.4 


-3.1 


- 


0.072i2 


0.03 


4.6 


-30 


— 


2.6 


-2.6 


- 


0.1822 


0.07 


4.6 


-43 


— 


2.6 


-1.7 


o.oe 23 


0.4122 


0.12 


4.5 


-31 


-185 


2.7 


-0.7 


0.28 23 


0.732 3 . 


0.18 


4.6 


-23 




2.9 


-0.9 


0.65 


1.10 


0.195 


4.7 


-22 




2.8 


-0.3 


0.74 


1.11 


0.215 


4.6 


-16 




2.7 


0.0 


0.91* 


1.14 


0.27 


4.6 


-11 


-155 


2.7 


0.1 




1.30 


0.38 


4.7 


-10 




2.8 


0.1 




1.22 


0.58 








3.4 


-0.1 




0.73 


0.68 


4.6 


3 




3.1 


-0.2 


1.06* 


1.14 


0.78 








3.6 


-1.9 


1.01* 


1.32 


1.0 


4.7 


-4 


-160 


3.6 


-1.8 


0.91* 


1.61 



this fit are given in Tab. |TTJ They are also plotted in 
Fig. [3] to show the overall evolution. It is clearly vis- 
ible that 9^' 2 4K (circles) is increasing monotonically 
with increasing x for x < 0.38 and even changes sign at 
x w 0.2. This leads to the important result that, with in- 
creasing x, the strength of the FM fluctuations increases 
and even dominate for 0.27 < x < 0.38 (9^' 2_4K > 0), 
although all samples of the series order antiferromag- 
netically at low temperatures. For the stoichiometric 
YbCo2Si2 compound the Weiss temperature is again neg- 
ative, 9* 2 ' 4K = (-1.8 ± 0.2) K. 

The second fit parameter of the low-temperature CW 
fit, ^^' 2_4K , reflects the evolution of the Kondo screen- 
ing. In Yb systems, increasing chemical pressure should 
lead to a decreasing Kondo temperature and, therefore, 
to a change in the dominating exchange interaction from 
Kondo to Ruderman-Kittel-Kasuya-Yosida (RKKY). In 
Fig. EJ the diamonds show the values of /i^Jj 2 ~ 4K for all 
x. The monotonic increase of /i^' 2 4K from 2.4 to 3.6 [Ib 
supports a decrease of the Kondo screening with increas- 
ing x, resulting in a larger effective moment. The final 
discussion about the interplay between Kondo effect and 
RKKY exchange interaction will be given in Sec. III. B. 
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FIG. 4: (Color online) Temperature dependence of the resis- 
tivity normalized to the room-temperature value for 0.18 < 
x < 1. The downward arrows mark the onset of coherence 
effect typical for Kondo lattice systems. With increasing x, 
Tmax shifts towards lower temperatures. The dashed arrow 
shows the CEF related hump for x > 0.58. In the inset, the 
low-temperature behavior of the resistivity for the Co con- 
centrations 0.18 < x < 0.27 is shown. Dash-dotted arrows 
indicate Tn 



C. Electrical resistivity 

In this section, we present the results of electrical resis- 
tivity measurements on single crystals with 0.18 < x < 1. 
The resistivity p(T) was measured between 0.35 < T < 
300 K using the four-point contact method and the ac- 
transport option of a Physical Property Measurement 
System (PPMS) from Quantum Design, equipped with 
a 3 He-insert. All samples were measured with the cur- 
rent perpendicular to the crystallographic c-direction 
and were normalized to their room-temperature value, 



p(300K), for comparison. In Fig. 21 this normalized re- 
sistivity is plotted on a logarithmic T-scale. For 0.18 < 
x < 0.68, a maximum is clearly visible (arrows in Fig. 0] 
are indicating the positions of the maximum, T max ). This 
results from the interplay between two different energy 
scales: (1) the CEF splitting and (2) the onset of the 
coherence effect related to the Kondo energy scale, Tk, 
as observed in various other Kondo lattices^ Since the 
susceptibility measurements show that the CEF does not 
change drastically with increasing chemical pressure, the 
shift of T max with x is unlikely to result from a change 
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of the CEF splitting. Furthermore, in previous pressure 
experiments on pure YbRh 2 Si2, T max also decreases with 
increasing pressure.— We, thus, conclude that the reduc- 
tion of T max for increasing x results from the decrease of 
Tk- A weak and broad hump in p(T) emerges around 
25 K (dashed arrow in Fig. 0]) for x > 0.58, which can 
be related to the CEF splitting. This clearly marks the 
detaching of the two energy scales, Tk and CEF effects, 
in p(T). For x = 1, T max is not any more distinguishable 
because the AFM order sets in before the coherence effect 
can develop. This leads to the result that at some point 
in the doping series the energy scales of the two interac- 
tions, Kondo and RKKY, must cross each other and the 
system changes from a more itinerant Kondo lattice to a 
magnetically localized system. A detailed discussion on 
this interplay will be given in Sec. III. B. 

Now, we will concentrate on the low-temperature be- 
havior around the AFM ordering. For all concentrations 
the onset of the antiferromagnetism at Tn can be de- 
duced from an anomaly in p(T) (see dash-dotted arrows 
in the inset of Fig. @| . Whereas for the concentrations 
0.27 < x < 0.58 and x = 1, the AFM ordering results 
in a drop of p{T), for 0.18 < x < 0.215 and x = 0.68 
we find an increase of p(T) just below Tn. This up- 
turn has already been observed for x = 0.07 (Refill) 
and 0.12 (Ref. 36 ) and probably emerges from a gap open- 
ing at the Fermi level below Tn in the new AFM Bril- 
louin zone. The second anomaly Tl for concentrations 
0.18 < x < 0.215 is again marked by a decreasing p(T). 
While this feature disappears with increasing Co concen- 
tration, it is recovered for x > 0.68. 

Analyzing the resistivity data in the vicinity of the 
magnetic transitions, another important fact can be seen 
in the inset of Fig. 2] Comparing the low-temperature 
behavior of /f(T) for x = 0.215 and 0.27, a pronounced 
change from a small increase to a distinct drop in p{T) 
at Tn can be noticed. This different temperature de- 
pendence is most likely due to a change of the magnetic 
structure between these two concentrations. Finally, we 
have to point out that for 0.18 < x < 0.215, the resistiv- 
ity is quasi- linear in T above the ordering temperature, 
as already shown for YbRli2Si2, 4 whereas for x = 0.27 
and x = 0.38 it is of a sublinear type (see inset of Fig. [4]) . 
A detailed analysis of the resistivity down to 20 mK will 
be published elsewhere^ 



D. Specific heat 

Now, we present the 4/ contribution to the specific 
heat, C 4 ? , for the single crystals with 0.18 < x < 1 at 
temperatures 0.35 < T < 5 K measured by the relaxation 
method in the PPMS in zero magnetic field. The lattice 
contribution has been subtracted from the total specific 
heat using the data for the non-magnetic reference com- 
pound LuRli2Si2. 37 Since the lattice contribution at 1 K 
amounts to less than 1% of the total measured specific 
heat, the synthesis and measurement of the appropriate 
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FIG. 5: (Color online) 4/ contribution to the specific heat for 
x = 0.18, 0.27, 0.58, and 0.78 evidencing different behaviors 
at the magnetic phase transitions in Yb(Rhi_ I Co I )2Si2. For 
x — 0.78, the corresponding relaxation curve at the first-order 
phase transition temperature, Tl = 1.01 K is plotted in the 
inset, with the arrow marking the onset of the latent-heat 
contribution. 



reference systems Lu(Rhi_, c Co a ;)2Si2 was not necessary. 

In Fig.[5j C 4/ /T vs. T for the concentrations x = 0.18, 
0.27, 0.58, and 0.78 is depicted. For x = 0.18, we detect 
two mean-held- type anomalies at Tn = 1.10 ± 0.03 K and 
T L = 0.65 ± 0.03 K. Below T L , the Sommerfeld coefficient 
7o = 1.3J/molK 2 is still very high. Raising the cobalt 
content to x = 0.27, the electronic specific heat changes 
significantly. For this amount of doping, C 4 ^/T shows 
only one A-type phase transition at Tn = 1.29 ± 0.03 K. 
The peak is much more pronounced and sharp, compared 
to the broad mean-field-typc anomalies for x = 0.18, 
which indicates, that the substitution of cobalt does not 
lead to significant disorder effects which would further 
smear the transition anomalies. The analysis of the crit- 
ical magnetic fluctuations around Tn reveals a conven- 
tional critical exponent for x = 0.27 and 0.38; 38 whereas, 
for x < 0.215 such an analysis could not be carried out 
due to the broader mean-field shape of the phase transi- 
tions. 

With increasing Co content the behavior of the spe- 
cific heat changes again. For x — 0.58, the AFM or- 
dering temperature Tn = 0.73 ± 0.03 K assumes a lower 
value, compared to the one for x = 0.38 (Tn = 1.22K), 
breaking the monotonic increase of Tn . This decrease of 
Tn between x = 0.38 and 0.58 will be discussed below. 
Further on, the absolute value of C 4 ^/T at Tn is en- 
hanced in comparison to the absolute values of the other 
concentrations and also the value at the lowest accessi- 
ble temperature of T = 0.35 K is considerably enlarged, 
which might be caused by another successive magnetic 
phase transition at lower temperatures. 

Two magnetic phase transitions are again visible for 
high concentrations of cobalt, e.g. x = 0.78, with 
T N = 1.32 ± 0.03K and T L = 1.01 ± 0.03K, as shown 
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FIG. 6: (Color online) 4/ contribution to the specific heat for 
all concentrations x > 0.18. The black and red lines mark the 
magnetic phase transitions Tn and Tl, respectively. 



in Fig. [5j The difference to the low-doped samples, e.g. 
x = 0.18, is that Tl shows clear signatures of a first- 
order phase transition. The values of C(T) across the 
sharp peak at Tl where taken continuously by calculat- 
ing the time derivative of the relaxation function plotted 
in the inset of Fig. The non-exponential temperature 
increase in time in the heating part is a direct evidence 
for a first-order phase transition. For x — 0.78, the upper 
transition at Tn = 1.32 K is a second-order mean-ficld- 
type AFM phase transition. It is, however, less sharp, 
probably owing to increasing disorder effects for large 
amounts of cobalt substitution. Furthermore, the Som- 
merfeld coefficient appears to be substantially smaller in 
comparison to the lower doped samples. 

The heat-capacity measurements of all doped sam- 
ples are depicted in a three-dimensional representation 
in Fig. [6] which displays the complexity of the magnetic 
behavior in this series. All extracted transition temper- 
atures are listed in Tab. HH Although several details re- 
garding the magnetic properties of Yb(Rhi„a; 00^)2812 
are not completely understood, a global picture may be 
drawn. The obtained results can be divided into three 
different regions depending on the amount of cobalt sub- 
stitution. (1) For 0.18 < x < 0.215 two phase transi- 
tion temperatures are clearly visible. With increasing 
Co concentrations, Tn is almost constant whereas Tl 
steadily increases towards higher values. In addition, the 
phase transition at Tl is changing its nature between 
x = 0.195 and x = 0.215 from second to first order. 
(2) Above x = 0.215, Tl merges the TN-Hne, since for 
0.27 < x < 0.58 only one transition temperature is ob- 
served. This phase transition is of second order, with 
Tn for x — 0.58 being reduced compared to the value 
at lower Co concentrations. (3) Above x = 0.58, a Tl- 
line again emerges from Tn . In contrast to those samples 
with lower Co content the transition at Tl is now always 
of first order and shifts towards lower temperatures with 



increasing Co doping. On the other hand, a second-order 
phase transition takes place at Tn which increases upon 
increasing substitution reaching a value of 1.61 ± 0.03 K 
for x = 1. Preliminary specific-heat measurements with 
applied magnetic field for selected concentrations show a 
shift of Tv and Tl to lower temperatures with increasing 
magnetic field confirming the antiferromagnetic nature of 
the phase transitions^ 



E. Photoemission spectroscopy 

The results of earlier measurements on YbCo2Si2 
performed by means of Mossbauer- and magnetic 
susceptibility^ techniques suggest that Yb reveals a 
purely trivalent state. On the other hand, for YbRli2Si2 
a slight deviation from trivalent Yb was found with a 
mean valence of 2. 9. 41 Hence, it is assumed that a grad- 
ual valence change can be followed when going from 
x = to 1. To this end, we applied X-ray photoe- 
mission spectroscopy (PES) to a series of single crystals 
with x = 0, 0.07, 0.12, 0.38, 0.68, and 1. PES addresses 
the electronic states under consideration directly and is 
thus an excellent probe to study the mean valence of a 
compound j 41 i 42 The experiments were performed at the 
SPring-8 facility using synchrotron radiation delivered by 
beamline BL25SU on single crystals cleaved immediately 
before the measurements, at 23 K and a pressure around 
2 ■ 10- 10 mbar. 

Fig. [7] presents the valence-band angle- integrated PES 
spectra taken at 450 eV photon energy for a series of 
Yb(Rhi_ x Co a ;)2Si2 single crystals. At this photon en- 
ergy PES can already probe deeper than the first few 
surface layers down to ~ 10 A,— ~— while the total en- 
ergy resolution of the instrument still remains reasonably 
high. The spectra exhibit two prominent regions, one be- 
tween ~ 1.3 eV binding energy (BE) and the Fermi level 
{Ep) emphasized by red, and a second one at 5 — 12 eV 
BE shaded light-gray. The spin-orbit split sharp doublet 
near the Fermi level represents the 4/ electron emission 
from divalent Yb atoms, while the more complex spectral 
structure at higher binding energies reflects the trivalent 
Yb stated The observed large energy shift between the 
divalent and trivalent Yb 4/ components is explained by 
the large Coulomb-correlation energy in the final state of 
photoemission. All spectra were taken at Si-terminated 
surface regions of the in situ cleaved samples. There- 
fore, contributions of purely divalent surface Yb atoms 
(broad doublet with components at 0.8 and 2.0 eV BE) 
are largely suppressed. 

The PES spectrum taken from pure YbRli2Si2 exhibits 
both divalent and trivalent 4/ components, typical for an 
intermediate-valent ground state of this material, with 
the possibility of strong hybridization between 4/ and 
valence states, which was indeed observed in our previ- 
ous experiments ! 41 ' 47 Upon gradual substitution of Rh 
atoms by Co, the intensity of the bulk divalent compo- 
nents close to Ep decreases considerably in comparison 
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FIG. 7: (Color online) Photoemission spectroscopy on cleaved 
single crystals for different x showing the gradual disappear- 
ance of the intensity of the bulk divalent components close to 
Ef with increasing doping. 



to the trivalent features and almost disappears for the 
pure YbCo2Si2 compound. Moreover, the trivalent mul- 
tiplet experiences a rigid shift towards the Fermi level 
by 0.45 eV. Reason for that might be a different chem- 
ical shift of the trivalent Yb component for YbRfi2Si2 
and YbCo2Si2. In fact, thermo-chemical calculations of 
the Yb 3+ energy shift for both compounds^ propose a 
shift to lower BE when going from the Rh to the Co 
compound. Moreover, the decrease of the divalent Yb 
component when going from x = to 1 shows that the 
respective state in the ground state shifts to higher ener- 
gies above Ep. Assuming that the Coulomb-correlation 
energy in the photoemission final state is comparable for 
both compounds the trivalent Yb component should shift 
towards Ep accordingly. Which of both effects is dom- 
inating here remains subject of further investigations. 
In any case, the experimental observations clearly re- 
veal the stabilization of the trivalent state of Yb ions 
in pure YbCo2Si2 which would characterize this material 
as showing a rather weak interaction between the 4/ and 
the valence states in comparison to YbRli2Si2. A thor- 
ough analysis of the mean Yb valency for the different 
concentrations based upon results of high-resolution res- 
onant inelastic X-ray scattering experiments, is presently 
underway. 49 



YbRh 2 Si 2 under pressure: ' * 
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2 4 6 8 10 12 14 16 18 20 
p (GPa) on YbRh 2 Si 2 

FIG. 8: (Color online) Magnetic phase diagram of the series 
Yb(Rhi_2,Coa:)2Si2 in comparison to the results from pressure 
experiments reported in the literature . 19 ' 25 ' 26 Clearly visible 
is the good agreement of the data for the low doped region. 
For high amounts of cobalt (x > 0.68) Tjv increases much 
faster under pressure than for Co substitution. 



III. DISCUSSION 

A. Magnetic phase diagram Yb(Rhi_ I Co I )2Si2 

The results obtained by magnetic, transport, ther- 
modynamic, and spectroscopic measurements lead to 
a consistent picture of the magnetic behavior of the 
Yb(Rh 1 _ x Co x )2Si2 single crystals. The x — T phase di- 
agram is presented in the inset of Fig. [5] For the sake 
of simplicity, the transition temperatures Tn and Tl are 
presented for the heat-capacity measurements only. For 
x < 0.12, the results from low-temperature susceptibil- 
ity data are taken from the literature^ First-order phase 
transitions at Tl for a; = 0.215 as well as for 0.68 < x < 1 
are marked by black circles. The transition temperatures 
for x > 0.18 obtained from resistivity and susceptibil- 
ity measurements are in very good agreement with these 
specific-heat results. Since the analogy between hydro- 
static and chemical pressure effects is one of the inter- 
esting aspects of this work, we compare in the main part 
of Fig. [H] our results for the alloy with the reported re- 
sults for pure YbRli2Si2 under pressure. By using the 
measured lattice parameters (see Tab. and the bulk 
modulus of YbRh 2 Si2 (189 GPa)^ we can convert the 
Co concentration into chemical pressure. 

We shall first address the phase diagram of the alloy 
series on its own, and then discuss the comparison with 
hydrostatic pressure results. The main features in the 
x — T magnetic phase diagram of the alloy are a steep 
increase of Tn(x) with x at low Co contents, a broad 
maximum near x — 0.3 with Tj^max = 1-30 K, a decrease 
towards a minimum near x = 0.58 with TN.min = 0.73 K, 
followed by a further increase towards TV = 1.61 K in 
pure YbCo2Si2- The transition at Tjv seems to be second 
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order in the whole concentration range. Further on, a sec- 
ond transition at Tl < Tv appears near x — 0.07, shifts 
continuously to higher temperatures with increasing x, 
changes from second-order to first-order type around 
x = 0.21, and finally merges with Tjv(x) near x = 0.27. 
A second transition at Tl < Tv, of first order type, is 
also observed at higher Co contents x > 0.58, beyond 
the minimum in Tv(x). Its critical temperature Tl is 
almost insensitive to the Co content. Our results indi- 
cate that in the whole concentration range the anomaly 
at Tjv corresponds to a transition towards an AFM state, 
despite the Weiss temperature becoming FM in the re- 
gion 0.2 < x < 0.5. We have no evidence for a change 
towards FM ordering. On the other hand, the anomalies 
at Tl for low or high Co contents are obviously con- 
nected with a change in the AFM structure, i.e., in the 
propagation vector and/or in the orientation of the Yb 
moments. This phase diagram is rather complex, with a 
pronounced non-monotonic behavior, indicating a com- 
petition between different interactions: e.g., Kondo ef- 
fect, different exchange paths, and CEF effects. 

It is therefore quite amazing that the main features 
of this complex phase diagram are nicely reproduced in 
the pressure experiments. The non-monotonous volume 
dependence of Tv is almost identical for Co substitution 
and hydrostatic pressure up to x — 0.58 or p = 8GPa. 
The appearance and initial increase of Tl is also iden- 
tical in pressure and Co substitution experiments, but 
the further evolution is different. In the alloy, Tl merges 
with Tv near x = 0.27, and we did not found any ev- 
idence for a second transition below Tn in the range 
0.27 < x < 0.58. In contrast, in the pressure experi- 
ments of Knebel et al£&, the second transition could not 
be resolved in the pressure range 2.3 < p < 4.6 GPa, 
but was again clearly observed in specific heat data at 
higher pressures, 4.6< p < 6.8 GPa. There, the transi- 
tion temperature Tl(p) first increases sharply and then 
stays slightly below Tv(p). 

The most pronounced differences arc however observed 
beyond x = 0.58 or p = 8 GPa. In diamond anvil 
cell experiments a sharp increase in Tn with p, up to 
Tn ~ 7.5 K is observed^ In contrast, in the alloy series 
we found a much weaker increase of Tjy, only up to 1.6 K 
in pure YbCo2Si2, and a further transition at Tl < Tn 
for x > 0.68. At least two effects can be proposed as 
origin of this clear difference: (1) The increase of c/a 
upon Co doping might result in a significant difference 
between the CEF under hydrostatic and chemical pres- 
sure at higher p. In this context one should note that 
the Mossbauer experiments in the diamond anvil cell in- 
dicate a switch of the direction of the ordered moment 
from the suspected in-plane direction at low p to along c 
for p > 9 GPa, while Mossbauer results on pure YbCo2Si2 
indicate the ordered Yb moment to be within the basal 
plane. 20 (2) For hydrostatic pressure only a small part of 
the volume reduction at p = 10 GPa can be attributed 
to the shrinking of the Yb atomic volume connected with 
the slight reduction of the Yb valence from « 2.9 at p = 




FIG. 9: (Color online) Evolution of the characteristic energy 
scales of Yb(Rhi_a;Co I )2Si2 illustrated by the parameters Tai , 
Tmax, and 70 evidencing the change of the dominant energy 
scale from Kondo to RKKY (dashed lines are guide to the 
eyes). The inset shows the increment of the entropy, AS, 
used to determine T±{. 

to pure trivalent at high pressures. The largest part of 
the volume reduction at higher pressures corresponds to 
a real compression of Rh, Si, and Yb atoms which should 
result in an increase of hybridization and thus an increase 
of the RKKY interactions. 



B. Kondo versus RKKY exchange interaction 

In the previous paragraph, the evolution of the mag- 
netic behavior of Yb(Rhi_ x Co x )2Si2 has been related to 
the change of the average unit-cell volume, V. One of the 
important mechanism for the understanding of the phase 
diagram is the competition between the local Kondo and 
the inter-site RKKY exchange interactions. This com- 
petition is directly connected to the strength of the hy- 
bridization between the localized 4/ states and the de- 
localized valence-band states, expressed by the exchange 
constant J which is a function of V. 

For the heavy- fermion metal YbRh2Si2, with a large 
Sommerfeld coefficient (70 = 1.7J/molK 2 ) this hy- 
bridization is strong and the Kondo effect dominates 
over the RKKY interaction. A Kondo temperature for 
the lowest-lying CEF Kramers doublet, Tk f» 25 K, has 
been deduced from the magnetic entropy 4^ For Yb com- 
pounds, the application of pressure leads to a stabiliza- 
tion of the magnetic, trivalent Yb state, which in a sim- 
ple Kondo-lattice model corresponds to a decrease of J. 
Due to the different analytical dependencies of the energy 
scales Tk and Trkky on the hybridization strength J, 
this decreasing J results in a much faster exponential de- 
crease of Tk compared with the quadratic decrease of the 
RKKY energy scale. 51 Hence, at some point the RKKY 
interaction becomes dominant, and the AFM ordering is 
stabilized. 



9 



The decrease of Tk and the change from a Kondo- to 
a RKKY-dominated behavior in Yb(Rh 1 _ x Co 2 ;)2Si2 is il- 
lustrated in Fig. [9] by plotting several parameters such as 
T 4 f, the Sommerfeld coefficient 70, and the position of 
the high-temperature maximum in p(T), T max , as func- 
tion of x. The parameter T 4 f (triangles in Fig- [9]) reflects 
the energy scale of all exchange interactions acting on the 
lowest CEF doublet and has been deduced from the mag- 
netic entropy, S, as T 4 f = 2 • T(S — 0.5i?ln2) (see arrows 
in the inset of Fig. Hat S = 0.5Mn2). For x < 0.38, T 4f 
is far above Tn, highlighting a dominating Kondo energy 
scale. With increasing amount of cobalt T 4 f strongly de- 
creases reflecting the decrease of Tk- However, the strong 
decrease of T 4 f stops for x > 0.4 and is replaced by an 
almost x-independent behavior. Since this change oc- 
curs where T 4 f is approaching Tv, it obviously reflects 
the crossover of the two energy scales Tk and Trkky, 
the later one becoming the dominant one for x > 0.6. 
This overall behavior is confirmed by the evolution of 
the coherence temperature T max (circles in Fig. [5]), de- 
duced from the maximum in p(T) (cf. Sec. Ill Cj ). At low 
Co concentrations, T max decreases strongly with increas- 
ing x, showing an almost perfect scaling with T 4 f. For 
x > 0.4, T max proceeds decreasing, in contrast to T 4 f, 
directly evidencing a further reduction of Tk- For x = 1, 
Tmax shifts below Tn, which implies Tk well below Tv. 
The decrease of Tk is nicely confirmed by the increas- 
ing effective moments at low temperatures, //^' 2_4K (see 
Tab, im and Fig. [3]), which reflects the diminished screen- 
ing of the local 4/ moments upon increasing Co content. 
The Sommerfeld coefficient 70 (squares in Fig. [5]) , which 
was deduced from a linear fit of C 4 ^/T versus T 2 below 
the phase transition temperatures, assuming a dominant 
contribution of AFM magnons to the specific heat, also 
decreases continuously with increasing Co content, from 
the huge value 70 = 1.7J/molK 2 in pure YbRh 2 Si2 to 
merely 70 = 0.13 ± 0.05 J/molK 2 in pure YbCo 2 Si 2 . The 
decrease of 70 is smoother than the decrease of Tk ■ From 
a very simple phenomenological approach one expects a 
decrease of 70 with Tk in the magnetic ordered region 
of the phase diagram of a Kondo lattice, but there is no 
strict theoretical prediction, and the relation is not uni- 
versal, since 70 depends not only on Tk, but also on the 
molecular field. 52 Further support for the decrease of the 
quasiparticle mass comes from the PES results for differ- 
ent x (see Fig. [7]) . The strong reduction of the density of 
states at the Fermi level due to a drop of the Yb 2+ peak 
at Ef, is directly visible. 

Finally, we discuss the 4/ contribution to the entropy 
which was obtained by integrating C 4 ^/T over T from 
the lowest measured temperature T = 0.35 K (cf. inset 
of Fig. [9J . Two important features can be deduced from 
the behavior of S(T) under increasing chemical pres- 
sure. First, since the entropy does not exceed Mn2 at 
T = 10 K throughout the whole series, only the lowest 
CEF doublet is involved in the formation of the ground 
state, the first excited CEF doublet remains energetically 
well separated. Second, the smooth and continuous in- 



crease of S with increasing amount of cobalt rules out 
an abrupt change from a low moment to a high moment 
state, but rather demonstrates a smooth stabilization of 
the AFM order. This behavior is evident from S(T) at 
Tn- For low x, S(Tys) is only a small fraction of i?ln2, 
e.g., S(x = 0.18, T = Tat) ~ 0.3i?ln2, but with increas- 
ing Co concentrations, S'(Tjv) rises continuously towards 
higher values, e.g., S(x = 0.78, T = T N ) ~ 0.7i?ln2. 



IV. CONCLUSION 

We have reported a detailed investigation of the evo- 
lution of magnetism in the alloy Yb(Rhi_ 2 ,Co a; )2Si2. We 
first described the crystal growth and the structural and 
chemical characterization, and then presented measure- 
ments of the magnetic susceptibility, electrical resistivity, 
specific heat, and photoemission for different concentra- 
tions covering the whole concentration range. The re- 
sults revealed a complex magnetic phase diagram, with 
an antiferromagnetic ground state at all compositions, 
but with a non-monotonic evolution of the Neel temper- 
atures Tjv and a second antiferromagnetic transition at 
Tl < Tjv for some concentrations. Tjv first increases 
with increasing x, from 72 mK in pure YbRh.2Si2 to a 
first maximum with Tjv = 1.3 K near x = 0.3, as ex- 
pected for an Yb system under chemical pressure. How- 
ever, for x > 0.3, instead of increasing further, Tv de- 
creases smoothly to Tv = 0.73 K near x = 0.58, and 
rises again smoothly to Tv = 1.6 K in pure YbCo2Si2. 
Some significant changes in the magnetic state seem to 
occur in the regions i«0.3 and x « 0.6, since they not 
only correspond to extrema in Tjv(x), but also to tricrit- 
ical points where a second transition line T^(x) merges 
with Tjv(x). The Tl lines extend from the two tricritical 
points at x w 0.27 and x « 0.58 to almost pure YbRh 2 Si2 
and pure YbCo2Si2, respectively. The anomalies at Tl 
at low or high Co contents are obviously connected with 
a change in the antiferromagnetic structure. While the 
transition at Tjv is always second-order type, the transi- 
tion at Tl is of first-order type in the whole region at high 
Co contents but switches from first-order to second-order 
type when moving away from the tricritical point at low 
Co contents. Despite its complexity, this phase diagram 
is very similar to the phase diagram under hydrostatic 
pressure when both are plotted as a function of the unit 
cell volume. For small reduction of the volume they are 
identical, while some differences appear at larger volume 
reduction, the most prominent being a strong increase of 
Tjv for p > 8 GPa contrasting a much smaller increase 
for large Co substitution. The strong similarity indicates 
that Co substitution basically corresponds to chemical 
pressure with no evidence for Co magnetism. Further- 
more, analysis of susceptibility, resistivity and specific 
heat data indicates a pronounced decrease of the Kondo 
temperature with increasing Co content. For x > 0.38, 
the Kondo temperature becomes comparable to Tn, re- 
sulting in a crossover from itinerant to localized Yb-4/ 
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states as well as a change of the dominant mechanism 
from the local Kondo to the inter-site RKKY interaction. 
This localization of the 4/ electrons is accompanied by 
drastic changes in the low-temperature properties studied 
by magnetization, specific heat, resistivity, and photoc- 
mission spectroscopy. In addition, the evolution of the 
susceptibility data in the series confirms a competition 
between ferromagnetic and antiferromagnetic exchange 
as first suggested from NMR data in pure YbRli2Si2. 18 
For 0.2 < x < 0.5 the Weiss temperature determined 
from Curie- Weiss plots at low T becomes ferromagnetic, 
although ordering stays antiferromagnetic. Therefore, 
the series Yb(Rhi_ a; Co 2: )2Si2 allows direct experimen- 
tal access to the two competing magnetic correlations 
in YbRli2Si2 at finite and zero wave vector and offers 
a unique opportunity to get a deeper understanding of 
the nearby unconventional quantum critical point. For 
example, tracing the magnetic propagation vector from 
YbCo2Si2 to YbRli2Si2 using neutron diffractometry may 
be one key experiment to clarify the magnetic structure 



of YbRh 2 Si 2 . 
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